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We apply the dynamical systems tools to study the linear dynamics of a self-interacting scalar field 
trapped in the braneworld, for a wide variety of self-interaction potentials. We focus on Randall- 
Sundrum (RS) and on Dvali-Gabadadze-Porrati (DGP) braneworld models exclusively. These mod- 
els are complementary to each other: while the RS brane produces ultra-violet (UV) corrections to 
general relativity, the DGP braneworld modifies Einstein's theory at large scales, i. e., produces 
infra-red (IR) modifications of general relativity. This study of the asymptotic properties of both 
braneworld models, shows - in the phase space - the way the dynamics of a scalar field trapped in 
the brane departs from standard general relativity behaviour. 
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I. INTRODUCTION 

Randall-Sundrum (RS) braneworld models have an ap- 
preciable impact on early universe cosmology, in particu- 
lar, on the inflationary paradigm. Actually, a distinctive 
feature of cosmology of a scalar field confined to a RS 
brane is that the expansion rate of the universe differs 
at high energy from that predicted by standard general 
relativity. This is due to a term quadratic in the energy 
density, that produces enhancing of the friction acting on 
the scalar field. This means that, in RS braneworld cos- 
mology, inflation is possible for a wider class of potentials 
than in standard cosmology [l[ . Even potentials that are 
not sufficiently flat from the point of view of the con- 
ventional inflationary paradigm can produce successful 
inflation. At sufficiently low energies (much less than the 
brane tension) , the standard cosmic behavior is recovered 
prior to primordial nucleosynthesis scale (T ~ 1 MeV) 
and a natural exit from inflation ensues as the field accel- 
erates down its potential Q. 1 Another interesting feature 
of this scenario is that the inflaton docs not necessarily 
need to decay; it may survive through the present epoch 
in the cosmic evolution. Therefore, it may also play the 
role of the quintessence field, which is a necessary ingre- 
dient to explain the current acceleration of the expansion 
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1 In this scenario, reheating arises naturally and radiation is 
created through gravitational particle production [J| and/or 
through curvaton reheating [J|. This last ingredient improves 
the brane "steep" inflationary picture Other mechanisms 

such as preheating, for instance, have also been explored [f|. 



of the universe. Such a unified theoretical framework for 
the description of both inflaton and quintessence with 
the help of just one single scalar field has been the target 
of some works (see for instance Refs. @, 0) B S (HI)- 

Another braneworld model that has received much at- 
tention in the last years, is the Dvali-Gabadadze-Porrati 
(DGP) model. It describes a brane with 4D world- 
volume, that is embedded into a flat 5D bulk, and al- 
lows for infrared (IR) /large scale modifications of grav- 
itational laws. A distinctive ingredient of the model is 
the induced Einstein-Hilbert action on the brane, that 
is responsible for the recovery of 4D Einstein gravity at 
moderate scales, even if the mechanism of this recovery 
is rather non-trivial [llj . The acceleration of the expan- 
sion at late times is explained here as a consequence of 
the leakage of gravity into the bulk at large (cosmologi- 
cal) scales, so it is just a 5D geometrical effect, unrelated 
to any kind of mysterious " dark energy" . As with many 
IR modifications of gravity, there are ghosts modes in 
the spectrum of the theory [12] • 2 Nevertheless, study- 
ing the dynamics of DGP models continues being a very 
attractive subject of research. It is due, in part, to the 
very simple geometrical explanation to the " dark energy 
problem" , and, in part, to the fact that it is one of a very 
few possible consistent IR modifications of gravity that 
might be ever found. 

Aim of this paper is to extend the study of references 
[l4l [Ti ] - the investigation of the dynamics of a self- 
interacting scalar field trapped on a DGP brane, and on 
a RS braneworld, respectivley - to include a wide vari- 
ety of self-interaction potentials beyond the constant and 



2 In fact there are ghosts only in one of the branches of the 
DGP model; the so called "self-accelerating" branch, or self- 
accelerating cosmological phase ll.'J . The Minkowski cosmologi- 
cal phase is free of ghosts. 
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exponential potentials. This goal will allow us to make 
conclusive arguments in favour of (or against) the claim 
made in (Iol | about the genericity of unification of the 
inflaton and of the quintessence in the Randall-Sundrum 
scenario. We expect a similar result regarding genericity 
of gravitational screening of the potential energy of the 
scalar field within the DGP brane context. 

As in references [HI, IH| here we make use of the dy- 
namical systems tools to retrieve useful information on 
the asymptotic properties of the models under study [l6| . 
In order to be able to analyze self-interaction potentials 
beyond the exponential one we will rely on a method 
proposed recently in the reference [l7| . 

The organization of the paper is as it follows. In sec- 
tion II we provide important details about the Randall- 
Sundrum model. These include the field equations, the 
phase space variables chosen, and the mathematical def- 
inition of the phase space itself. The same features, this 
time for the Dvali-Gabadadze-Porrati braneworld model, 
are given in section III. The results of the study of the 
corresponding critical points and their stability proper- 
ties are shown in section IV. The section V is aimed at 
the physical discussion of the above results, while the 
conclusions are given in section VI. Through the paper 
we use natural units (87rG = &Tr/m Pl = h= c = 1). 



II. THE RANDALL-SUNDRUM MODEL 

We will be concerned here with the dynamics of a self- 
interacting scalar field with an arbitrary self-interaction 
potential, that is trapped in a Randall-Sundrum brane of 
type 2 (RS2). The field equations - using the Friedmann- 
Robertson- Walker (FRW) metric - are the following: 
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where A is the brane tension, 7 is the barotropic index of 
the background fluid, pr = Pc/> + Pm and V is the scalar 
field self-interaction potential. 

Following [l5} we introduce the following dimensionless 
phase space variables in order to build an autonomous 
system out of the above system of cosmological equations: 
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where the comma denotes derivative with respect to the 
new time variable (r = In a). Notice that, after the above 
choice of variables one can realize that 
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This means that the four-dimensional (4D)/low-energy 
limit of the Randall-Sundrum cosmological equations - 
corresponding to the formal limit A — > 00 - can be as- 
sociated with the value z = 1. The high-energy limit 
A — > 0, on the contrary, corresponds to z — > 0. The 
critical points associated with z = 0, if any, have to be 
analyzed carefully. Actually, in connection with the clas- 
sical character of the underlying theory of gravity, the 
physical meaning of these points in phase space has to 
be taken with caution due to the high energies associ- 
ated with them. 

As long as one considers just constant and exponen- 
tial self-interaction potentials {d^V = and d^V = 
const respectively), the equations ((5][7|) form a closed 
autonomous system of ordinary differential equations. 
However, if one wants to go further to consider a wider 
class of self-interaction potentials beyond the exponential 
one, the system of ordinary differential equations ([5][7]) is 
not a closed system of equations any more, since, in gen- 
eral, 8$ V is a function of the scalar field 0. A way out 
of this difficulty can be based on the method developed 

in E3- 

In order to be able to consider arbitrary self-interaction 
potentials one needs to consider one more variable s, that 
is related with the derivative of the self-interaction po- 
tential through s = —d^V/V = — d^lnV. Hence, an 
extra equation 



-V6xs 2 (T - 1) 



(9) 



V 



V3H' 



Pt 

m 2 ' 



After this choice of phase space variables we can write 
the following autonomous system of ordinary differential 
equations (ODE): 



has to be added to the above autonomous system of equa- 
(4) tions. The quantity T = Vd 2 V/ (d,pV) 2 in equation ((U is, 
in general, a function of tfi. The idea behind the method 
in [13] is that T can be written as a function of the vari- 
able s G 5R + , and, perhaps, of several constant param- 
eters. Indeed, for a wide class of potentials the above 
requirement: T = T(s), is fulfilled, see Tabled] 
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As in [13 we introduce a new function f(s) — T(s) — 1 
so that equation ^ can be written in the more compact 
form: 



s' = -V6xs 2 f(s). 



(10) 



Equations (f5T l7ll0p form a four-dimensional closed au- 
tonomous system of ordinary differential equations, that 
can be safely studied with the help of the standard dy- 
namical systems tools [la ]. 

The phase space for the autonomous dynamical system 
driven by the evolution equations (|5l l7ll0p can be defined 
as it follows: 



* = {(x,y, z) : < x 2 
-1 < x < 1,0 < y,0 < z < 1} x { s e 3t + }. (11) 



III. 



THE DVALI-GABADADZE-PORRATI 
MODEL 



In this section we will focus our attention in a 
braneworld model where a self-interacting scalar field is 
trapped on a DGP brane. In the flat FRW metric, the 
field equations are the following: 



(12) 
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model with constant and exponential potentials, exclu- 
sively. The study of other potentials was not considered. 

As in the former section, following a method devel- 
oped in [It]], here we extend the analysis of the three- 
dimensional autonomous system of reference [l4j , to four- 
dimensions, through the addition of the extra- variable 
s = —d^V/V defined in the former section. This will 
permit us to consider a wider class of self-interaction po- 
tentials beyond the exponential one. In consequence, to 
the system of equations of [l4j, we add (p~0|) . so that we 
are left with the following autonomous closed system of 
ordinary differential equations: 
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where, as before, the comma denotes derivative with re- 
spect to the time variable r = / Hdt, and f(s) — T(s) — 1 
(T^VdlV/^V) 2 ). 

After the above choice of phase space variables, the 
equation (fT4")) can be put into the following form: 



where p m is the energy density of the background 
barotropic fluid (7 is its barotropic index), V its self- 
interaction potential, (f> is the scalar field trapped in the 
DGP brane, and 
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H z ± — H, 



(14) 



where, as customary, r c is the crossover scale inherent in 
the DGP brane model. There are two possible branches 
of the DGP model corresponding to the two possible 
choices of the signs in (fT4"|) : "+" is for the Minkowski 
cosmological phase of DGP model - that is free of ghosts 
-, while "-" is for the self-accelerating solution. 

Following the reference [3] we define the phase space 
dimensionless variables: 



V 
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The corresponding autonomous system of equations in 
the variables x, y and z, defined above, was used in [14j to 
study the asymptotic properties of the DGP-quintessence 
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For the Minkowski phase, since < H < 00 (we con- 
sider just non-contracting universes), then 1 < z < 00. 
The case — 00 < z < — 1 corresponds to the time rever- 
sal of the later situation. For the self-accelerating phase, 
— 00 < z 2 < 1, but since we want real valued z only, then 
< z 2 < 1. As before, the case — 1 < z < repre- 
sents time reversal of the case < z < 1 that will be 
investigated here. Both branches share the common sub- 
set (x,y,z — 1), which corresponds to the formal limit 
r c — > 00 (see equation (|14p ). i.e., this represents just the 
standard four-dimensional behaviour typical of Einstein- 
Hilbert theory coupled to a self-interacting scalar field. 

The phase space for the autonomous system f| 164119(1 . 
for the "+" branch can be defined as: 



3 In fact, fitting SN observations requires H > r^ 1 in order to 
achieve late time acceleration (see, for instance, reference [l2| 
and references therein). This means that z has to be real- valued. 
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TABLE I: Explicit form of the function V for several quintessential potentials. 

r(s) Potential Reference 

1 + ± - s£ V = Vo sinh- Q (A0) [18] 

s 2 + 2ms 2 +8mA + s^/ s 2 +8mA -tr- -ct- exp[A</> 2 ] r.. q-i 

1 - h + # V = % [cosh(A^) - If [20] 

_ K[Ka p+s (a+m v = Vq [exp(QK0) + ex p(/3 K 0)] [21] 

-| 7 = Vo exp(-\<j>) + A [22] 



TABLE II: Properties of the critical points for the autonomous system (|5I7I10|1 . 
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while, for the self-accelerating "-" phase, it is given by 
the non-compact region: 



9- = {(x,y,z): 0<x 2 +y 2 <l, 

z G]0,1]} x {s e 5R+}. (22) 

Notice that the points belonging in the set (x, y, 0) can 
not be included since, in this case (z = => Q = 0), 
the variables x and y are undefined. The self-accelerating 
solution H — l/r c (Q_ =0 z = 0) has been studied 
in [3] • In that reference the analysis of the critical points 
of the quintessence model under study was based on a 
concrete form of the self-interaction potential. Here, as 
in the former section, we use the approach proposed in 
[TtJ to investigate the critical points of the dynamical 
system for arbitrary functions /(s), so that, in principle, 
we are able to study arbitrary self- interaction potentials. 



IV. CRITICAL POINTS OF THE 
AUTONOMOUS DYNAMICAL SYSTEM 

In this section we will analyse in detail the crit- 
ical points of the autonomous systems corresponding 
to both Randall- Sundrum and Dvali-Gabadadze-Porrari 
braneworld models, as well as their stability properties. 



A. The Randall-Sundrum braneworld 

The critical points of system (|5l l7ll0|) are summarized 
in table|n] The eigenvalues of the corresponding jacobian 
matrices are shown in table IIIII In both cases s* is the 
value which makes the function /(s) vanish, i. e., /(s*) = 
0. In the same way we have chosen 

as 

As we see from tableslllland lIII| the point P\ exist in all 
cases regardless of the form of the self-interaction poten- 
tial (arbitrary s). Points P2 — Pi are associated with po- 
tentials whose first </>-derivative vanishes at some/several 
point /points (this case includes the constant potential 
whose (^-derivatives at any order vanish everywhere). 
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TABLE III: Eigenvalues for the critical points in table HI A = ^{2 - 7)(247 2 - s^(9 7 - 2)) 
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Worth noticing that the existence of points P^ 1 — P7 de- 
pends on the concrete form of the potentials (recall that 
the s»-s depend on the functional form of f(s)). From 
the table of eigenvalues, notice, besides, that there are 
four non- hyperbolic critical points/sets of critical points 
(at least one of the eigenvalues is vanishing): Pi, P2, P3, 
and P4. 

We recall that four-dimensional effects are associated 
with points belonging in the plane (x,y,z = 1). For 
points with z 7^ 1, five-dimensional effects affect the 
dynamics of the universe. There is only a set of criti- 
cal points with z / 1 (represented by P2 in table [IT} : 
(x, y, z, s) = (0, y, y 2 , 0). For points in this set, since 
x = 0, z = y 2 , while pr = V, then the Friedmann equa- 
tion can be written in the form: 



3P 2 = W 1 + — J . (23) 

For values of the potential much larger than the brane 
tension V A => Prs = V/V&X, so that the early 
time/high energy expansion rate in the Randall-Sundrum 
model (Hrs) gets enhanced with respect to the general 
relativity rate: 



H RS 
H G r 




This is the way brane effects fuel early inflation in the 
RS model. The fact that this is a critical point in the 
phase space of the RS model means that, helping infla- 
tion to happen is a generic feature of Randall-Sundrum 
braneworld models. For a further discussion about this 
see [H. 

The rest of the critical points of the dynamical sys- 
tem ([51 171101) lie on the plane (x,y, 1), so that, only four- 
dimensional behavior can be associated with them. Non- 



hyperbolic critical points in the set Pi = (0, 0, 1, s) corre- 
spond to the matter-dominated solution and, as it is seen 
from table IIII1 these are always saddle points in phase 
space. Critical points P2 — P4 have been exhaustively 
studied in [14j . 

Points are saddle critical points. These correspond 
to the solution dominated by the kinetic energy of the 
scalar field (f2^ = 1). This result differs from the one 
in standard four-dimensional theory, where the kinetic 
energy-dominated solution can be a past attractor (an 
unstable source point) for trajectories in the phase space. 

The non-hyperbolic point P3 has a two-dimensional 
stable subspace, which corresponds to a late-time attrac- 
tor solution (3P 2 = V). 

There are other two critical points that can be associ- 
ated with late-time attractor solutions: P§ and P7. For 
values s 2 < 37 (s*d/ > 0), the scalar field-dominated 
solution (point P§) is the future attractor of the au- 
tonomous system (|51 l7ll0p . The scaling solution (point 
P7) is the late-time attractor whenever it exists. 

Due to the fact that there are several non-hyperbolic 
critical points which can not be consistently studied 
with the help of the present linear analisys, we choose 
several concrete examples - corresponding to different 
quintessence potentials of cosmological interest - in order 
to illustrate, in the phase space, the dynamical behavior 
of the corresponding RS model, including the neighbour- 
hood of these non-hyperbolic critical points. 



1. V = V {smh(\(p)]- a 

This potential was studied for the first time in [24j . 
where it was shown to be a new cosmological tracker 
solution for quintessence. For this potential the function 
/(s) = 1/a — a\ 2 /s 2 , while 
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s* = ±aX, df = 



2aX 2 



(25) 



The scalar field-dominated solution (set of points Pq 
in tables [Til and IIII[) is a late-time attractor whenever 
A 2 < 37/a 2 , and a > 0. It is consistent with accelerated 
expansion if A 2 < 2/a 2 , and a > 0. Since s* = ±aA, P7 
is a stable node: 
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Therefore, whenever the condition 37/a 2 < A 2 < 6/a 2 
(a > 0) holds true, the matter-scaling solution domi- 
nates the late-time evolution of the universe. In fig[T] 
this behaviour is illustrated for an arbitrary set of initial 
conditions. For a > 0, and A 2 > the matter-scaling 
solution is a stable spiral. It is the late-time attractor in 
this case. 
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This behavior is clearly shown in fig. [JJ It is seen, in 
particular, that the trajectories in phase space emerge 
from the point S = (x,y,z) = (0,0,0) - the empty 
Misner-RS universe - meaning that this is the past at- 
tractor of the Randall- Sundrum cosmological model. We 
want to notice that the points with z = have been re- 
moved from the phase space ^ since, in general, at z = 
the autonomous system of equations (|5 l7ll0p blows up 
due to our choice of phase space variables. For that rea- 
son the point S does not appear in table UJl See [l5[ for 
further discussion. 



FIG. 1: Trajectories in phase space (x,y,z) for different sets 
initial conditions for the potential V = Vb[sinh(A0)] -Q1 . The 
free parameters have been chosen to be (0,7, A): (3, 1,0.5) - 
upper panel, (3,1,0.65) - central panel, and (3,1,5) - lower 
panel. For the first parameter selection the late-time attractor 
is the scalar field-dominated solution (point P 6 in table 
For the remaining parameter selections the late-time attractor 
is the matter-scaling solution (point P 7 in table Hill- In the 
last case this point is a stable spiral. 



2. V = F o [cosh(A0) - If 

This potential was proposed in [25j in order to describe 
both quintessence and a new form of dark matter called 
frustrated cold dark matter, due to its ability to frustrate 
gravitational clustering at small scales. Accelerated ex- 
pansion is obtained for < a < i. For this potential 
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This potential (first studied in [26|) drives the evolu- 
tion of the universe to transit from a scaling attractor 
into a de Sitter-like attractor. Following the above ex- 
plained methodology we have: 



-a/3, df = --j. 



(28) 



(29) 



The condition s*df > is satisfied whenever /3 > 0, 
so that the possible late-time attractor solutions are the 
following: 



If - J §J < a < y Ti 
P-j is a stable node 



!> the matter-scaling solution 



• If a > 



stable spiral. 



the matter-scaling solution P7 is a 



As in the above examples, from figd] we see that the 
past attractor in the phase space is the empty, Misner-RS 
universe. 



B. The Dvali-Gabadadze-Porrati braneworld 



FIG. 2: Phase trajectories (x,y,z) for given initial data, for 
the potentials V = Vb[cosh(A</>) — l] a (upper panel), and V = 
(, ) + e -°a.»)ij (lower panel). The past attractor in the phase 
space is the empty, Misner-RS universe (point S = (x, y, z) — 
(0,0,0) that does not appear in table [TTJ) . In the upper panel 
we set a = 0.4, A = 5 and 7 = 1 so that the scalar field- 
dominated solution (point Pe) is the late-time attractor. In 
the lower panel the matter-scaling solution Pi is a stable spiral 
(a = -2.6, f3 = 2 and 7 = 1). 



while 



s* = ±a\, df = -. (27) 

si 

In order for the critical points Pq and P7 to be late-time 
attractors, it is necessary that the condition s„df > be 
fulfilled. For the present potential this condition can be 
written in the following form: 




This constraint is fulfilled only for negative a < 0. There- 
fore Pq and P7 are both saddle points. The critical point 
P3 is the late-time attractor (see fig. [5]). 



We study here the stability of the critical points of the 
autonomous system (|16H19[) . As before, we will rely here 
on the assumption that the function f(s) has zero(s) at 
given value(s) s*: /(s*) = 0. 

The critical points of the system (|16H19[) are shown 
in table IIVI while the eigenvalues of the corresponding 
Jacobian matrices are shown in the table El We have not 
included the points with z = —1 in our analysis, since the 
equations (|16H19j) are invariant under the change of sign 
z — > —z. 

The critical points P± , P 2 ± , P3 are associated with the 
stationary points of the potential - extrema and saddle 
stationary points - and, in general, with potentials whose 
(^-derivative vanishes (including the constant potential). 
The existence of the remaining critical points depends, in 
general, on the concrete functional form of the potentials, 
since, as stated before, the value s* is determined by the 
form of f(s). 

The points Pi,P 2 ± ,P 3 ,P 4 are non-hyperbolic critical 
points (one of the eigenvalues of the corresponding Ja- 
cobian matrixes vanishes). In this case the only thing 
we can state with certainty, on the basis of the straight- 
forward analisys of the autonomous system of equations 
(|16til9[) is that, depending on the phase considered - the 
Minkowski phase or the self-accelerating one -, and on 
the initial conditions, trajectories in phase space origi- 
nating in one of the repeller points (P^ 1 ), will inevitably 



8 



TABLE IV: Properties of the critical points for the autonomous system (|16H19I 



Pi 


X 


y 


2 


S 


Existence 




O»0 


q 


Pi 








1 





both branches 





undefined 


37-2 
2 


p ± 


±1 





1 







1 


1 


2 


P 3 





i 


2 







1 


-1 


-1 


Fa 


A 

u 


U 


1 


s* 


always 


U 


undefined 


37-2 
2 


pt 


±1 





1 


s* 


always 


1 


1 


2 


Pe 


s * 

v/TT 




1 


s* 


s* < 6 


1 


"2-3 
3 


2 


P< 


3 _7_ 

2 s„ 


/ 3 7 (2-7) 
V 2^ 


1 


s* 


always 


37 

7T 


-1 + 7 


-1 + 22 
1 ~ 2 



TABLE V: Eigenvalues for the critical points in table HVTA = sj{2 - 7)[24 7 2 - s|(9 7 _ 2)] ). 



P l 


Ai 


A 2 


A 3 


A 4 


Pi 





1(7 "2) 


37/2 


3 7 /2 


p ± 


3 


3 





3(2 - 7) 


Ps 


-3 








-37 


Pa 





1(7-2) 


37/2 


37/2 


P? 


3 


+^6s» df 


3(2-7) 


3 T |s* 


Pi, 


2 


-*2 4f 


si - 37 


2 


Pi 


-3s* df 


^-[,.( 7 -2)-A] 


J-[ S .( 7 _2) + A] 


37/2 



approach one or several of the above non-hyperbolic crit- 
ical points. 

Points Pi and Pa represent the matter-dominated so- 
lution, while and P 5 ± are associated with the solu- 
tion dominated by the kinetic energy of the scalar field 
(the stiff- matter solution). 4 They are linked always with 
decelerated expansion (q — 2). The de Sitter-DGP - 
accelerated - solution corresponds to the critical point 
P3. The point Pq, which exists whenever si < 6, is as- 
sociated with the scalar field-dominated phase, while P7 
represents, in the phase space, the matter-scaling solu- 
tion. 



Maximum we can say about the non-hyperbolic points 
is that they have attached an unstable subspace that is 
spanned by the eigenvectors: 



In fact, the points Pi and P 2 are particular cases of P4 and of 



''1 




1 

voy 



v 2 



1 





v 3 



(1\ 




voy 



P 5 , respectively, when s* 



0. 



The point P 5 is always a past attractor in phase space, 
while P+ is a saddle in ^. This is the classical result 
within general relativity with a minimally-coupled (self- 
interacting) scalar field. 

The scalar field-dominated solution (point P@) an d the 
matter-scaling solution (critical point P7), are always 
saddle points in the phase space. This result has to be 
contrasted with the classical result within general rela- 
tivity with a minimally-coupled scalar field, where, de- 
pending on the values of the constant parameters, the 
above mentioned solutions can be late-time attractors in 
the phase space. 

The above mentioned results are illustrated in the fig- 
ures [3] and [5]) for several quintessential potentials. 

For the potential V = V [sirih(\<t>)]~ a [HI (upper 
panel in the figure [3J = ±aA and df — +2/(a 2 A), 
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FIG. 3: Trajectories in phase space (x, y, z) for different 
sets of initial conditions for the self-accelerating phase of 
the DGP model The upper panel is for the potential 

V — Vo[sinh(A<^)] -t *, the panel at the center is for the poten- 
tial V = Vo exp(A0 2 ) / cj> a , while the panel at the bottom is for 
the potential V = Vb[cosh(A0) - If . 



FIG 
El 



4: Flux in time (x(t),y(t),t) for the potentials in figure 



s* = ±V~ SmX. The bottom panel is for the potential 



V = Vo[cosh(A</«) - If [20]. 5 In this case one have that, 



while for V = Vq exp(A0 2 )/</> m [l9[ (a potential origi- 
nated in supergravity models - middle panel in fig. [3]) 



5 The author of reference [20ll found that, for small values of p, 



10 



at f(s) = 0, s* = ±Xp. 

Phase trajectories in (figure [3]) originate from the 
source critical points P^ 1 , corresponding to the standard 
4D kinetic energy dominated (stiff-matter) solution and 
(asymptotically) approach to the point (0,1,0,0) that 
has been removed from the phase space since phase space 
variables x and y blow up at the phase plane (x, y, 0, s). 
The dynamics in the neigbourhood of this point has to be 
investigated in terms of different phase space variables. 

Phase trajectories in (figure [5J originate from the 
4D stiff-matter solution (unstable node P 5 ± in tab. fV|) 
and end up at the inflationary points (0, 1, Zoi, 0) £ P^, 
where the different zoi-s are associated with the different 
initial conditions. Otherwise, points in P^ 1 are seen as 
attractor points by the different phase space " observers" , 
moving along different phase trajectories that originate 
at P+ or P~. 

V. RESULTS AND DISCUSSION 

A. RS model 

The main results of section IV. A can be summarized 
as follows: 

• The matter-dominated solution (point Pi) is a non- 
hyperbolic critical point independent of the func- 
tional form of the self-interaction potential. It can 
be, at most, a saddle. 

• For si < 37 (s*df > 0), the scalar field-dominated 
solution (point P§) is the future atractor. Other- 
wise, Pq is a saddle point. As seen from table |TT] 
this critical point can be associated with acceler- 
ated expansion whenever si < 2-f (s*df > 0). 

• For values si > 37 the matter-scaling solution 
(point P7) is a the late-time attractor. For 37 < 
si < 6 it is a stable node, while, for si > 6 it is a 
stable spiral. This solution is always a decelerating 
one. 

• The kinetic energy-dominated/stiff fluid solution 
(point Pjf) is always a saddle critical point in the 
phase space. 

The non-hyperbolic critical point P2 (in fact, a set of 
critical points), represents the slow-roll Friedmann equa- 
tion relating the Hubble expansion parameter with the 
potential of the inflaton field, modified by the presence 
of the RS brane (see equation (|23|) ). 




FIG. 5: Trajectories in phase space (x, y, z) for different 
sets of initial conditions for the Minkowski phase of the 
DGP model *]/+. The upper panel is for the potential V = 
Vo[sirih(A</>)]~ Q , the panel at the center is for the potential 
V — Vo exp(\(f> 2 ) /(f> a , while the panel at the bottom is for the 
potential V = Vo[cosh(A<^) - l] p . 



(p < 1/2), the scalar field dominates the mass density in the 
universe at late times, leading to accelerated expansion. This 
potential might serve as a good candidate for quintessence. The 
matter-scaling is approximately constant during the prolonged 
epoch. 



In general, the dynamical behavior of the Randall- 
Sundrum model differs from the standard behavior 
within four-dimensional Einstein-Hilbert gravity coupled 
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to a self-interacting scalar field, only at early times (high 
energy regime). Actually, the empty (Misner-RS) uni- 
verse is always the past attractor in the phase space of the 
Randall-Sundrum cosmological model This result is 
to be contrasted with the standard four-dimensional re- 
sult where the kinetic energy-dominated solution is the 
past attractor [23|, [26| . Within the present scenario the 
latter solution (critical points P4 and P^ 1 in tables HT1 and 
IIIID is always a saddle point. 

The late-time cosmological dynamics, on the contrary, 
is not affected by the RS brane effects in any essential 
way. 

B. DGP model 

^From the analysis in section IV. B, the following im- 
portant results can be summarized: 

• Points Pi — P4 in table IIVI are non- hyperbolic crit- 
ical points so that, only the behaviour of the phase 
space trajectories may uncover the main properties 
of the dynamical system in their neighbourhood. 

• The kinetic energy-dominated solution can be ei- 
ther a past attractor (point P 5 ~) or a saddle point 
in the phase space as in [23| - see figure[21 where the 
time evolution of the dynamical system is shown - 
since, at early times, the DGP brane effects can 
be safely ignored so that the standard cosmologi- 
cal dynamics is not modified. Recall that the DGP 
brane effects produce infra-red modifications to the 
laws of gravity. 

• The scalar field-dominated solution (critical point 
Pq in Tab. IIV[) . as well as the matter-scaling so- 
lution (point P7), represent always saddle points 
in the phase space, contrary to the classical re- 
sult within four-dimensional general relativity plus 
a minimally-coupled scalar field. 

Apart from the exponential potential, there are a large 
number of potentials that can produce the matter-scaling 
solution (critical point P7). This result is expected since, 
in the 4-dimensional limit, when standard Friedmann be- 
havior is recovered, we are left with the case studied in 
the reference (23j . where the matter-scaling solution was 
identified as a critical point in phase space. 

Nevertheless, the DGP brane effects indeed modify the 
late-time cosmological dynamics through changing the 



stability of the corresponding (late-time) critical points. 
Actually, in the present case the matter-scaling solution 
(critical point P7), as well as the scalar field-dominated 
phase, are always saddle critical points. This result has to 
be confronted with the classical general relativity result 
where the above mentioned solutions can be a late-time 
attractors. 

VI. CONCLUSIONS 

In the present paper a thorough study of the phase 
space of both, the Randall-Sundrum and the Dvali- 
Gabadadze-Porrati braneworlds - with a self-interacting 
scalar field trapped on the brane -, has been undertaken. 
A wide class of self-interaction potentials for which the 
quantity T = VdiV/(d$V) 2 can be written as a function 
of the variable s = —d^V/V, are included in this study. 

It has been demonstrated, in particular, that the 
empty Misner-RS universe is always the past attractor 
in the phase space of the Randall-Sundrum cosmolog- 
ical model. The RS brane effects modify the early- 
time dynamics, so that, additionally, the kinetic energy- 
dominated solution - the past attractor within gen- 
eral relativity plus a slef-interacting (minimally-coupled) 
scalar field - is always a saddle critical point. The critical 
points that can be associated with late-time behaviour, 
as well as their stability properties, are not modified by 
the RS brane effects. 

A detailed study of the dynamics of the DGP brane 
(with a self- interacting scalar field trapped on it), re- 
veals that the critical points in phase space, coincide 
with the ones found in standard (four-dimensional) gen- 
eral relativity. An additional critical point that can be 
associated with five-dimensional behaviour, can be found 
only for the constant self-interaction potential. Neverthe- 
less, even if, in general, there are no critical points that 
could be associated with genuine higher-dimensional ef- 
fects, DGP brane effects indeed play a role: they modify 
the stability properties of the critical points associated 
with late-time cosmological dynamics. 

The above results have been clearly illustrated with 
the help of phase space pictures generated by several po- 
tentials of cosmological interest. 
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of the present research. 
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